Sources and holes in a one-dimensional traveling wave convection experiment 
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We study dynamical behavior of local structures, such as sources and holes, in traveling wave 
patterns in a very long (2 m) heated wire convection experiment. The sources undergo a transition 
from stable coherent behavior to erratic behavior when the driving parameter e is decreased. This 
transition, as well as the scaling of the average source width in the erratic regime, are both qualita- 
tively and quantitatively in accord with earlier theoretical predictions. We also present new results 
for the holes sent out by the erratic sources. 

PACS numbers: 47.20.Bp, 47.20.Dr, 47.54.+r, 07.05.Kf 



Traveling wave systems play an exceptional role within 
the field of pattern formation. If the transition to pat- 
terns is supercritical (forward), the dynamics close to 
threshold should be amenable to a description by the 
Complex Ginzburg-Landau (CGL) amplitude equation 
0. Theory and experiments are difficult to compare, 
however, for the following two reasons: (i) both the CGL 
model and experimentally observed traveling wave pat- 
terns exhibit an astonishing variety of ordered, disordered 
and chaotic dynamics, which can be difficult to character- 
ize or compare, (ii) The dynamics depends, in general, 
strongly on non- universal coefficients H, but the 
values of these coefficients are difficult to determine in 
experiments [1,|EI3- 

The study of local structures, such as sources, fronts 
and holes which play an important role in traveling wave 
systems ^ iSi, ^ Jd^ , provide a promising 

route to compare theory and experiment as they partially 
circumvent these difficulties: their nontrivial behavior of- 
ten depends only on a subset of the coefficients 11] and 
is, in addition, relatively easy to characterize experimen- 
tally HSilll. 

In this paper we present a successful example of this 
approach in a heated wire convection experiment (Fig.^l. 
This system forms left and right traveling waves that sup- 
press each-other; typical states consist of patches of left 
and right traveling waves, separated by sources (which 
send out waves) and sinks (which have two incoming 
waves) Earlier theoretical work ^3|> which was 

based on the amplitude equations (|1I2|I below, predicts 
that, essentially due to the transition from an absolute 
to a convective instability |l5j |. sources tend to display a 
diverging width when the driving parameter e is lowered 
beyond a critical value [Eq. ©]. More recently it was 
predicted that just before these stationary sources would 
diverge, they become unstable and give way to fluctuat- 
ing sources of finite average width which display highly 
non-trivial dynamics 0. 

We indeed observe this nontrivial change in source be- 



0.7 m 




2 mm 



■2m- 



2. 5 cm 

FIG. 1: Schematic side-view and cross section of the heated 
wire experiment. A thin (0.2 mm diameter) wire (1) is 
stretched beneath the free surface of a fluid (depth h = 2 mm) . 
When it is heated by sending an electrical current through it, 
surface waves are excited. The slope of the waves is measured 
by reflecting a laser beam off the surface onto a position sensi- 
tive detector (PSD) (3). The laser and the PSD are mounted 
on a cart (2) that rides on precision steel rods (4). 



havior when the driving (heating of wire) is decreased; 
not only the measured transition value, but also the qual- 
itative behavior of sources is in accord with the predic- 
tions 4. 14] . All properties necessary to compare theory 
and experiment are measured in a set of independent ex- 
periments. The fluctuating sources send out holes, and 
we show that these display behavior very similar to that 
predicted for homoclinic holes Q. 



I. EXPERIMENTAL SETUP 

Our experiment consists of a 2 m long heated wire of 
diameter of 0.2 mm and resistivity of 50 fl m^^ that 
is placed under the free surface of the fluid at a depth 
h ~ 2 mm (see Fig. The wire is stretched to the 
breaking limit and its maximum sagging is 0.1 mm. The 
heat Q dissipated in the wire drives the system; through 
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a combination of gravity- and surface tension induced 
convection, surface waves emerge at Q = Qc that travel 
along the wire |0| . The sides of our cell are made of brass 
and contain copper tubes through which cooling water of 
21.0±0.1°C is circulated. In order to guarantee a clean, 
uncontaminated free surface, we use a low- viscosity, low- 
surface-tension silicon oil 

A sensitive linear measurement of the surface slope 
along the cell is obtained by recording the reflection 
of a laser beam off the fluid surface onto a position 
sensitive device. Both laser and position detector are 
mounted on a computer-controlled cart which travels on 
precision stainless steel rods. This allows us to mea- 
sure surface wave amplitudes as small as 0.5 /im. The 
signals of the scanning device are wave frequency and 
Hilbert transformed to yield the complex valued field 
A{x^t) = |A|exp(i0). From this the local wave num- 
ber is computed as q{x,t) = d(j>{x,t)/dx. To improve 
the signal to noise ratio, running averages over a time 
interval of 10 s are performed. 

Vince and Dubois already demonstrated that the 
primary bifurcation at Q = Qc is supercritical and ex- 
plored the phase diagram as a function of Q and wire 
depth h. For e ^ 0.15 the amplitude exhibits the scaling 
1^1 ~ £^/^. This is expected near a supercritical bifur- 
cation, and it also sets the range of applicability of the 
amplitude description. 



II. AMPLITUDE EQUATIONS 

For systems with counter-propagating waves, the 
appropriate amplitude equations are the coupled one- 
dimensional CGL equations 0: 
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Here Ar and Al are the amplitude of the right and left 
moving waves, sq is the linear group velocity, and ci,C2 
and C3 measure the linear and nonlinear dispersion. The 
experimentally accessible control parameter e measures 
the distance from threshold. The coefficients to,^o and 
go give the scales of time, space and amplitude. To model 
our experiment, where left and right traveling waves sup- 
press each-other, g2 should be larger than go 0. 



A. Scaling 

Sources show complicated behavior within the ampli- 
tude equations (HEJ Hill. For 
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FIG. 2: Numerical results for the behavior of sources in 
the coupled amplitude equation Q, lEJ- (a) Inverse aver- 
age source width as a function of e, for the coupled CGL 
equations with sq = 1.5, ci = —1.7, C2 = 0, ca = 0.5, go = 1 and 
g2 — 2. The coefficients so and ci were chosen to be similar 
to those measured in the experiment; also 172 > go in the ex- 
periment. The values of C2 and C3 where chosen such that the 
plane waves remain stable; their precise value does not play 
a significant role then. Note the crossover near Sca ~ 0.14. 
(b) Space-time plot of the local wave number of a fluctuat- 
ing source for £ = 0.11 < ~ 0.14, illustrating fluctuations 
of the width. In the black region the amplitude has fallen 
below 10% of the saturated value; the light and dark curves 
correspond to hole-like wave number packets send out by the 
source. 



sources are coherent structures with a well-defined shape, 
while for e < sources start to fluctuate and their aver- 
age width scales as oc (see Fig. 12)). The quantity Sca 
in Eq. (O is simply the value of e where the transition 
from absolute to convective instability of the A = state 
occurs ^3,01; its relevance can be understood as follows. 
Consider the dynamics of a single front in the left-moving 
wave amplitude Al only, for which Al{x ^ 1) ^ 0. The 
propagation velocity of this front is given by a compe- 
tition between the linear group-velocity, which tends to 
convect any structure to the left with velocity sq, and the 
propagation of the front into the A = state w ith ji the 
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comoving frame, velocity v* = 2^0 V^U + 
the front velocity in the lab frame is thus v* 
ing a source as a pair of fronts in Al (on the left) and Ar 
(on the right), it is clear that these fronts move together 
when e > Sca, but move apart when e < Sca] the change 
in direction of front propagation precisely corresponds in 
the transition from absolute to convective instability. 

Numerical simulations of Eq. ^ were done in order to 
see whether the experimentally observed source behavior 
described below could be understood on basis of the am- 
plitude description. Such simulations f3| have revealed 
that sources do not simply move apart and diverge when 
the instability of the A = state becomes convective; for 
e = e^" ^ Sea, when the sources have become very wide, 
they start to fluctuate. For smaller e, the average source 
width scales as (see Fig. 2). The mechanism respon- 
sible for the sources staying at a finite but large average 
width is not completely understood and may depend on 
the noise strength. In the low noise limit, the "tip" re- 
gions of the two fronts sense the other mode which leads 
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FIG. 3: Determination of the coefficients of the CGLE. (a) 
Timescale r determined from the exponential growth of the 
amplitude in quench experiments for various value of the heat- 
ing power Q. This time scales as r = ro/e, with tq =16(1) s. 
(b) Correlation length ^o- Full line: histogram of squared 
modulus \A\^ vs q which is measured from a modulated wave 
field at £ = 0.10. Dashed line: fit of \A\^ = 1 - Co/e(<3 - lof, 
with 5o = 2.7(6) X 10~^ m. \A\ is normalized so that |^| = 1 
corresponds to waves with wave number qo- (c) Front veloc- 
ity. Shown is the modulus |A(a;,t)]. The a;— extent of the scan 
is 0.682 m, the total time is 5242 s. At t = the power is 
quenched from Q = to Q = (1 -|- e)Qc, with e = 0.051. The 
white lines outline two fronts. Since e < e"° here, both fronts 
propagate in the same direction (here to the right). 



to the formation of phaseslips there. The resulting per- 
turbations are then advected by the group-velocity and 
amplified by the linear growth-rate, resulting in a jittery 
motion of the front. For larger noise strength, convective 
amplification of noise may compete with this mechanism. 

These phenomena are illustrated in Fig. |21 which has 
been calculated for parameter values which are in the 
range of the experimental ones, but we emphasize that 
the predicted source instability is generic and insensitive 
to the precise parameter values. 



III. MEASUREMENTS 
A. Front and group velocities 

Now that we have discussed the theoretical predictions 
for sources, we return to our experiment. For a compar- 



ison of the source behavior with theory, we need to de- 
termine the transition from convective to absolute insta- 
bility, which requires measurements of the group velocity 
and the front velocity as function of e. 

The group velocity sq was determined from the prop- 
agation of deliberate perturbations of the surface. We 
found that it has the same sign as the phase velocity and 
that is shows only a weak g— dependence, so we associate 
the measured value 2.1(1) x 10^** ms^^ with the linear 
group velocity sq. 

Fronts where made by quenching the heating power Q 
to a finite value Q ^ {1 + e)Qc at t — 0. After a short 
while, waves invade the unstable surface in the form of 
fronts. The boundaries of these fronts travel with SQztvf 
respectively, where Vf scales with e as u/o^/e. Fig. 
shows the evolution of the amplitude of the waves for 
e = 0.051; this value appears to be below Sca because the 
velocity of the fronts has the same sign as the group ve- 
locity. The results of several experiments, both at e < Sca 
and e > Sca yields that w/o = 5.4(5) x 10^^ m/s. Com- 
paring this to the value obtained for sq, we immediately 
find that Sca ~ 0.15(5). A alternative estimate of Sca 
was made from observing at which e the slowest moving 
edge of a front has zero velocity, which led to a value of 

Eca = 0.10(2). 



B. Measurements of the coefficients 

In principle, a confrontation of theory and experiment 
can also be performed by measuring the characteristic 
time (ro) and length (^o) scales, the linear group velocity 
So and linear dispersion coefficient ci, since from these 
the transition from convective to absolute instabilities 
also follows (see Eq. 3). Note that starting from the full 
hydrothermal equations, these coefficients can in princi- 
ple be obtained from a systematic amplitude expansion 
[l|. At present, we can only obtain ci via measurements 
of the front-velocity which leads to a consistency check 
(see below). The length and time scales are relevant for 
comparing spacetime diagrams of experiment and theory 
and are measured independently. 

The characteristic time is determined from measure- 
ments in which the growth of the amplitude is followed 
when, after a sufficient long transient in which a plane 
wave is established, e is increased from e ~ 0.017 to 
larger values. The initial growth of |^| is exponential 
~ exp(i/r), and repeating this experiment for various 
values of e yields the data presented in Fig. 13 Using 
that T scales as r = ro/e, we obtain that tq =16(1) s. 

The length scale was measured from weakly mod- 
ulated waves in the single- wave domains; according to 
Eq. (Q) with Al = 0, these are related by A{q)'^go — 
(e — Co9^)- Plotting these values we obtain Fig.|2Ib), in 
which we recognize the quadratic behavior of \A\ as a 
function of q. The measured differed slightly (but not 
systematically) from run to run and from a series of such 
measurements and fits we find for the correlation length 
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FIG. 5: (a) Wave number profile of a fiole emitted from the 
unstable source of Fig.UJj. Dashed line: typical wave number 
profile. To help reading wave number off the vertical axis, 
the plot has been sectioned, (b) Scatter-plot of the mini- 
mum of the modulus versus the extreme (in x) of the wave 
number along each of the holes shown in Fig. 0]d. Both com- 
pression (large q) and dilation (small q) holes belong to a 
one-parameter family. 



FIG. 4: (a) Dependence of the width of a source on the re- 
duced control parameter e. Dots: mean width {w)~^ , dashed 
line: fit of (w) ~ . (b) Dependence of the rms fluctua- 
tion a = V< > — < ui >2 of the width of a source on e. 
Notice that the source becomes unstable for e > Sca- Note 
that in (a)-(b), length scales have been non-dimensionalized 
by the characteristic scale fo. (c) Space-time diagram of the 
wave number field q{x,t) of an unstable source, e — 0.11; the 
extent of the a;— axis is 158 ^o, the total time is 660 tq. 

^0 — (2.7 ± 0.6) X 10^^ m, which only close to threshold 
becomes similar to the basic wavelength of the traveling 
waves. 

Taking these time and length scales and the measured 
front velocity Vfo used before, we find then that (l-|-c^) w 
2.6, from which it follows that ci = ±1.3(4). A weak 
consistency check is that [l + cf) should be larger than 
one; independent measurements of ci would lead to a 
stronger consistency check. 



01)- Fig. 0] shows that the behavior of (w) as a func- 
tion of e in the experiments shows the same qualitative 
features as the numerical simulations of Fig. |21 For de- 
creasing e the width appears to diverge, but at e « 0.15 
there is a cross-over to a {w) ^ e^^ behavior. Below 
this crossover value, the sources fluctuate strongly and 
the standard deviation of the width rapidly increases 
[l7| . In a cyclic fashion, these sources first grow and 
spawn outward-spreading wave fronts, leaving an inter- 
val of near-zero wave amplitude behind in the source core. 
Here phase slips occur, which make the fronts jump back; 
the resulting phase twists are carried away by hole struc- 
tures which travels roughly with the group velocity (the 
light and dark lines). In our numerical simulations of 
of the coupled CGL equations (Fig. |2Jl exactly the same 
hole structures are observed. The crossover value for £*° 
is consistent with the transition value Sca that we deter- 
mined before. 



IV. COMPARISON BETWEEN EXPERIMENT 
AND THEORY 

A. Sources 

Now that all relevant parameters of the amplitude 
equations are approximately known, we turn to the be- 
havior of sources in our experiment. The dependence of 
the width w{t) of sources on the control parameter e was 
measured in long experimental runs in which a source 
was located at large heating power e « 0.3 and then fol- 
lowed at progressively smaller values of e At each e, 
the source was observed for several hours by scanning the 
fiuid surface, while keeping the experimental conditions 
constant. 

From the width w(ti) at discrete scan times ti we com- 
puted the mean (w) (as well as the standard deviation a 



B. Holes 

The structures sent out by the erratic sources display 
a dip in the amplitude \A\ and are therefore referred to 
as holes. It is well known that holes play an important 
role in the dynamics of traveling wave systems, and that 
different types can be distinguished by whether the wave 
numbers of their two adjacent waves are similar or sub- 
stantially different 0, IS S [s", Tl] . From the measured 
wave number profile in Fig. [SJl it can be seen that the 
wave numbers at the back and front side of the hole are 
very similar. We therefore associate these holes with so- 
called homoclinic holes Q. In addition, they display the 
following typical homoclinic hole behavior (see Fig. 2b 
and 4b): they do not send out waves and occur quite 
close together, they can evolve to defects and their prop- 
agation velocity (which in lowest order is given by sq) 
depends on the value of the extremum of q. In the lo- 
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cal wave number plot of Fig. [3 dilation holes (the dark 
lines) have a larger velocity than compression holes (light 
lines), just as in the experimental plot Fig. 0];. In fact, 
the correlation between the type of wave number modula- 
tion and the velocity of these coherent structures depends 
on the sign of ci, which was selected accordingly for the 
numerical simulations. 

Since homoclinic holes are dynamically unstable, their 
local profiles slowly evolve along a one-parameter fam- 
ily; on a scatter plot of the values of the minimum oi \A\ 
versus the correspondingextremum of q, these values col- 
lapse on a single curve The holes in our experiment 
precisely show this behavior: The extrema of |v4| and q 
rapidly evolve toward a parabolically shaped curve, and 
stay there during their further evolution (Fig. Eb). We 
only observed these holes in our experiment for at most 
a few characteristic times — too short to see clear signs 
of the weak instability predicted from the CGL equation. 

V. DISCUSSION AND OUTLOOK 

Our experiments raise a number of suggestions for fur- 
ther work, (i) The width where sources start to fluctuate 



is larger in the theory (0(100 ^o)) than in experiments 
(0(20 ^o))j while the fluctuations appear stronger for ex- 
perimental sources. Experimental noise or non-adiabatic 
effects which perturb the fronts may play a role here. 
(ii) Earlier experiments 12] have shown that for differ- 
ent heights of the wire, qualitatively different behavior 
occurs. Systematic measurements of the coefficients as a 
function of height may turn the heated wire experiment 
into a CGLE-machine with tunable coefficients. (Hi) 
Longer observations and more controlled generation of 
holes may shed more light on their relation to the homo- 
clinic holes predicted by theory, and may show the highly 
characteristic divergence of lifetime as a function of ini- 
tial condition (iv) Sinks show non-adiabatic phase- 
matching and in fact posses completely anti-symmetric 
profiles QJ^; there is no clear theoretical understanding 
of this. 
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